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I. INTRODUCTION 



Integrability vs Quantum Thermalization 

Jonas Larson 
Department of Physics, Stockholm University, AlbaNova University Center, Se-106 91 Stockholm, Swede^ 

Non-integrability is often taken as a prerequisite for quantum thermalization. Still, a generally 
accepted definition of quantum integrability is lacking. With the basis in the driven Rabi model we 
discuss these issues. The only symmetry of the model is the total energy and it would be classified 
as non-integrable according to the most commonly used definitions. Despite this fact, a thorough 
analysis conjectures that the system will not thermalize. Thus, our findings suggests: (i) care should 
be paid when linking non-integrability with thermalization, or {ii) standardly used definitions for 
quantum integrability are unsatisfactory for certain cases. 
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CsJ ■ The concept of integrability is well defined in classical systems jll, '2'| . Integrability here means that the number of 

degrees of freedom is smaller than the number of independent constants of motion. Constants of motions in classical 
systems are characterized by vanishing Poisson brackets, and independence between them by mutually vanishing 
' r^ ' ' Poisson brackets. Classical integrability implies that the solutions are periodic and live on a torus of constant energy t3!| 
rS |, in phase space. Translating the above definition to quantum Hamiltonian systems directly leads to complications and 
. K ■ there is no accepted definition of integrability in quantum systems [J|. As an example, the number of degrees of 
C^ , freedom for a quantum systems is generally taken as the dimension of the Hilbert space, and in particular the Hilbert 
2 ' space dimension can be finite. Classically, the degrees of freedom are necessarily continuous variables and there seem 
to be a contradiction in having a well defined quantum-classical correspondence, i.e. finite size Hilbert space systems 
do not have a proper classical limit. The spin, for example, is a pure quantum property. 

With the development of techniques in isolating and controlling quantum systems [^-Q , questions regarding quan- 
tum integrability have gained renewed interest. Of special interest is out of equilibrium dynamics in closed quantum 
^ ' systems [7]. Cold atom systems are especially practical for in situ measurements of quantum many-body systems 
and they thereby also provide a handle to study pure quantum evolution [8]. As a result, long standing questions 
in quantum statistical mechanics can now be addressed in an experimentally controlled way. Of particular interest 
^fs is the long time evolution and whether an interacting quantum system equilibrates and if so what characterizes the 
• I relaxed state [9] . A common believe is that for a non-integrable quantum system the state thermalizes, by which we 
f^ [ mean that expectations of any local observable A can be evaluated from a micro-canonical state pMC- This conjec- 
^f) ' ture has been supported in several numerical studies of various models JOj, JJj]. However, it seems strange to coin 
such an assumption based on a concept that still today lacks a proper definition. Moreover, it has been numerically 
demonstrated that using standard definitions for quantum integrability one can find models that are non-integrable 
• y—t ' and still do not thermalize J 12| . 
p\^ , Quantum thermalization has become deeply connected to interacting many-body systems [7|. It is important 

JH ■ to understand, however, that there is nothing in the theory that relies on having a quantum many-body system, 
- - -' i.e. a system possessing many degrees of freedom. It is rather properties of the eigenstates and the spectrum that 
determine the fate of the state 9]. Indeed, quantum thermalization has been demonstrated in systems whose classical 
counterparts possess only two degrees of freedom [13| . In the works of ref. [l3|, a common feature is instead that 
the corresponding classical models are chaotic [ij|. A question thereby rises: Is classical chaos a common feature of 
systems that quantum thermalize? Naturally, this cannot be a general condition since, as argued above, some quantum 
systems do not have a well defined classical limit. For example, the work [12] considers a disorder Heisenberg spin- 
1/2 chain for which a proper classical limit does not exist. To spur the discussion about quantum integrability and 
thermalization, in this work we consider a purely quantum model (i.e. lacking a classical counterpart) that does not 
obey standard criteria for integrability and still do not show any signatures of thermalization. More precisely we 
analyze out of equilibrium long term evolution in the driven Rabi model (RM) which describes interaction between a 
spin-1/2 system and a boson mode. After a general discussion about integrability and thermalization we investigate 
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statistical properties of various local expectation values. Despite the fact that the driven RM appears non-integrable 
according to standard definitions, all our results indicate absence of quantum thermalization. 

II. INTEGRABILITY AND QUANTUM THERMALIZATION 
A. Quantum Integrability 

Already mentioned in the introduction, integrability in classical systems has a clear meaning. An A^-dimensional 
Hamiltonian system 7J(p,q) is said to be integrable if: (i) there exist N single- valued constants of motion /„, i.e. 
{In,H} — 0, where { , } denotes the Poisson bracket, (m) the constants of motions /„ are functionally independent, 
and (Hi) the constants of motion /„ are in involution meaning {/„,/„'} = 0, V n,n'. For an integrable system, 
the solutions (p(t),q(i)) are periodic and evolve on {N — l)-dimensional tori in phase space. For such constrained 
evolution, the solutions do only explore a small part of the phase space. When the integrability condition is (slightly) 
lifted, the tori start to deform in accordance with KAM-theory (Kolmogorov-Arnold-Moser) [2|. The solutions are (in 
general [15|) no longer periodic and cover a larger part of phase space. This describes the transition from regular to 
chaotic motion in classical systems. 

Trying to define integrability for a Hamiltonian quantum systems H is far from trivial [J, Il6l Il7| . There have 
been numerous different attempts to give a meaningful and consistent definition. We summarize some of them in the 
following list. 

1. Traditional I. The far most commonly used definition for quantum integrability is obtained from translating 
the classical definition into a quantum language. Thus, functions /„ are replaced by operators /„ and Poisson 
brackets by commutators { , } — > z[ , ]/h. It is easy to reject such a definition by noticing that the projectors 
Pa = \'>Pa){ipa\, with \ipa) a (non-dcgenerate) eigenstate of the Hamiltonian define constants of motion and 
mutually commute. Thus, it is possible to find a set of constants of motion such that any quantum system 
appears integrable. 

2. Traditional II. The problem with the above definition led to the notion of relevant and irrelevant constants of 
motion 16, ISj. The relevant constants of motion are those which can be associated with a classical counterpart. 
This again have some flaws since inequivalent quantum constants of motion can share the same classical limit jl6| , 
and not all quantum systems do have a classical limit to start with. 

3. Scattering. A quantum system is integrable if its scattering is non-diffractive [17| . This applies only to continuous 
models and it relies on properties of the asymptotic scattered states. Thinking in terms of a scattering problem, 
if the outgoing solution contains "diffrective contributions" the system is non-integrable. 

4. Bethe solution. A quantum system is integrable if it can be solved with the Bethe ansatz. 

5. Poissonian level statistics. A quantum system is integrable if its energy level statistics is Poissonian \\^ . 
Following ref. [l9|, this definition relies on semi-classical arguments and to systems with continuous degrees of 
freedom. 

6. Level crossings. A quantum system is integrable if it shows level crossings. This definition is related to the 
previous one since avoided crossings are characteristic for systems showing level repulsion^ i.e. the energy level 
statistics follows a Wigner-Dyson distribution |14l | . Note that the definition does not say anything about avoided 
crossings. 

7. Solvability. A quantum system is integrable if it is exactly solvable. 

It can be argued that the defining properties of (iii), (v) and (vi) are rather consequences of non- integrability than 
defining it. The usefulness of definitions (iv) and (vii) may be discussed (for obvious reasons). We should mention 
that the list above is not complete, there exist further definitions not included here. 

B. Quantum Thermalization 

In recent years we have seen an increased interest in dynamics of closed quantum systems 7] . An open question 
with a very long history concerns equilibration of such states [2^ . A central topic in this field has been to understand 
local relaxation to a thermal state of a quantum many-body state |9l-ll2l. l2ll | . To gain deeper insight in the mechanism 



driving quantum thermalization, several concepts have been introduced, for example: eigenstate thermalization hy- 
pothesis (ETH) [9J] , quantum central limit theorems [22| , system-bath entanglement [23i | and eigenstate randomization 
hypothesis J24l |. 

Especially the ETH has been thoroughly studied. To explain the idea of ETH, let us express some time evolved 
state in eigenstates Itpa) of the Hamiltonian, 

|vl'(t)) = ^ae-^"*/''|^.). (1) 

1/ 

The expectation of some observable A reads 

(i(i)) = J2 \Cc\^A,, + ^ C:C^e'(^--^-)*/''^,^, (2) 

where A^^ = (■(/'i,|A|^^). If the state equilibrates, the long time expectation {A)^'^ should attain the time averaged 
value 



(i)^T ^ hm 1 rdt{A{t)) = J2 |C.I'^- 



(3) 




This expectation is obtained when the long time state is diagonal in the eigenvalue basis /5lt — ^2^ \Cv\^\tjJiy){tpiy\- 
For the situations of interest for us, the probabilities jC^P are only non-zero in some energy window AE around 
E = ($|i/|^). The number of populated states \ijju) in the sum ^ can be estimated with the inverse partition 
ratio [25| 

/ \ ~^ 

(4) 

Clearly, 77^ ^ 1 in order to expect equilibration. For a microcanonical distribution, /5mc — N{E,d)^^ J^jes l''!'-/) ii^il 
where S (> AE) is again an energy window around E and N{E, 5) being the number of states within (5, the expectation 
become 

The state thermalizes if {A)^^ = {A)^'^ up to corrections of the order 0(77!^). Now, the ETH says that for a state 
that thermalizes, Aaa varies little within the energy interval AE. We directly see that if A^a is more or less constant 
in the interval of interest the expectation {A)^"^ approximates the one you get from a microcanonical distribution. 
Thus, ETH predicts that any functions {A) has a weak .©-dependence on the scale AE. 

The ETH does not say whether a system will thermalize or not, it is rather a property of a system that thermal- 
izes [26|. Without deeper reflection, it is often assumed that any non-integrable system will thermalize. From the 
discussion in the previous subsection it is clear that there is a great ambiguity in such an assumption, simply because 
there is no generally accepted definition of quantum integrability. The problem might be circumvented for systems 
with a well defined classical limit, and it has indeed been found that several systems with a chaotic corresponding 
classical counterpart do thermalize. While numerical experience indicates such a fact, there is no strict proof that this 
is true in general. The situation is more complicated when the system of interest does not allow for a simple classical 
limit. 

Of course, the discussion above on the ETH is fully general, i.e. there are no assumptions on number of degrees 
of freedom nor on existence of a classical limit. As the name suggests, it relies on the properties of the eigenvectors. 
In the next section we will study a particular model which should not be classified as integrable according to the 
definitions above, and still we find no indications of quantum thermalization. 

III. GENERALIZED RABI MODEL 

A. Driven Rabi Model 

The RM 2t\ has a long history in quantum optics and especially in cavity quantum electrodynamics (QED) [28|. 
Despite its simplicity, a spin-1/2 system coupled to a single boson mode, the physics is extremely rich. In most 



experiments to date, both in cavity and circuit QED, the rotating wave approximation (RWA) is weU justified and 
the RM can then be solved exactly in terms of the Jaynes-Cummings model [29l] . Within the RWA, the number of 
excitations N — a^a + az/2 (a^ and a are the creation and annihilation operators for the boson mode and az is the 
Pauli ^-matrix acting on the spin) is preserved which implies a continuous U{1) symmetry, i.e. [t/^, Hjc] = with 
U^ = e*^"^. More recently, an alternative RWA was considered in order to derive an analytically solvable model with 
a larger validity regime compared to the regular RWA [30|. Also this time, the applied approximation results in the 
same U{1) symmetry. Relaxing the RWA means that the U{1) symmetry is broken down to a discrete Z2 parity 
symmetry {[U^, Hr\ = 0) characteristic for the RM [3ll]. In the spirit of the previous section, it is not clear whether a 
discrete symmetry results in integrability of the RM. Furthermore, while the boson mode has a well defined classical 
limit the spin does not and one cannot thereby define integrability from any classical limit. 

The search for a solution of the RM has a long history [32l- A breakthrough came in 2011 when D. Braak claimed to 
have solved the RM [33.] . In particular, the spectrum can be divided into a regular and an exceptional part. The regular 
part is given by zeros of a transcendental function. The exceptional solutions have a simple analytical expression but, 
on the other hand, t hey only exist for certain system parameters. More recently, A. Moroz remarked that the RM is 
not exactly solvable |34| . but rather an example of a quantum model that is quasi- exactly solvable [35:]. Thus there is 
a debate whether the RM is in fact integrable/solvable or not. 

We may break the Z2 symmetry by considering the driven RM, 
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9\ 



Xa^r- 



(6) 



Here, we have introduced dimensionless parameters by letting the energy hil of a single boson set a characteristic 
energy scale, uj is the energy separation of the two spin states ]1) and ]2), a^ is the Pauli x- matrix, g is the spin-boson 
coupling, and finally A is the drive amplitude. By letting A = we regain the Rabi Hamiltonian Hfj. The drive term 
breaks the parity symmetry since llT^axU^^ = —(Jx- Note that the drive of the spin can be removed by unitarily 

transform the Hamiltonian with the displacement U = e ^9 ^ -^ . In return, the transformed Hamiltonian contains 
a drive of the boson mode, i.e. [a^ + a) X/g. Thus, driving of the spin or the boson mode is unitarily equivalent and 
here we consider the first option. Judging from refs. [33] and [Sl], it seems that also the driven Rabi model dU is of 
the quasi-exactly solvable type. This fact may naturally be of importance in terms of thermalization. 

Let us return to the definition of quantum integrability in the previous subsection and check whether the driven 
Rabi model fulfills any of them. 
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Figure 1: The two adiabatic potentials V^{x) for the parameters lo = 1, g = 10 and A = 2. The inset shows a zoom of the 
avoided crossing. 



1. Traditional I. As already pointed out, this definition is pointless since one can always find a set of constants of 
motion such that any system would be considered integrable. 

2. Traditional II. With the driving, the Z2 symmetry is broken and the only relevant constant of motion is the 
energy. In this respect, the driven Rabi model should not be classified as integrable. Of course, we have 
a problem here since the spin does not have a natural classical limit. We may, however, perform a semi- 
classical approximation in which the boson mode is treated at a mean-field level, while the spin is still kept 



as a quantum entity. Thus, we make a coherent state ansatz for the boson field where the bosonic operators 
are replaced with their corresponding coherent amplitudes, a — ?> a and a) — ^ a*. In doing this we neglect 
any quantum correlations between the spin and the boson mode. As a result [3a|, a generic spin state can 



be written |6) 



i±z_ 

2 



^e'^* 



(tan(A0) = {dy) / {(Jx))- By introducing quadratures x andp according to a 
we can write a "classical" Hamiltonian 



where Z is the inversion {{(Tz) = Z) and A^ the relative phase 

{x + ip)/^/2 and a — {x~ip)/^/2, 



n X oj 



The semi-classical equations of motion now become 



x=p, 



gxV2 + x\ v^l - Z^cos{A^ 



gV2VT^^cos(A^), 



(7) 



Z = gxV^Vl - Z^sm{A^), A^ = ^ - (gV2x + x] cos(A0) 



Z 



(8) 



vr^z2 



Putting A = we obtain the classical equations of motion for the Dicke model which have been demonstrated 
to be chaotic [13, [33] ■ The corresponding semi-classical equations of motion for the RM were also analyzed in 
ref. [3^ with clear signatures of chaos. See also ref. [S^ which studies similar semi-classical equations of motion. 
We have solved the equations of motion (|8]) numerically and studied different Poincare sections [40l| . For large 
enough couplings g they all show well developed chaos. In this respect, the RM should not be considered 
integrable. 
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Figure 2: Level statistics of the driven Rabi model for the dimensionless parameters uj — 1, g — 10 and A = 2. Energies 
< E < 250 have been considered. 



3. Scattering. Since the spectrum of H is discrete, the idea of non-diffractive scattering does not apply to our 
system. 

4. Bethe solution. The Bethe ansatz is typically applied to quantum many-body problems with continuous degrees 
of freedom. Hence, we cannot apply such approaches to the RM. 

5. Poissonian level statistics. Level statistics explores the distribution P{S) - the number of energies with certain 
nearby energy gaps Sn = En+i — En. Typical for systems showing regular dynamics is that the level statistic of 
the spectrum follows a Poissonian distribution P{S) = e~ . Chaotic systems, on the other hand, show instead a 
level-repulsion effect and normally a Wigner-Dyson distribution P{S) = {S'n:/2)e^^''"/^ [14|. The level statistics 
of the RM has been studied in the past [41]. Despite the similarity to the Dicke model, their statistics are very 
different. While the Dicke model shows clear level- repulsion in the chaotic regime [37|, the level statistics of the 
RM is neither of Poissonian nor Wigner-Dyson shape. This was also pointed out by D. Braak in [33] where he 
noticed that the energies are rather equally spaced throughout. Many of t he p roperties of the spectrum can 
be understood within the Born-Oppenheinier approximation (BOA) 28, 4^|43|. In the BOA we decouple the 
internal degrees of freedom of H by diagonalizing the spin part of eq. ^ . The two resulting adiabatic potential 



VtA^)^^±\-r+ V29X + X) . (9) 



surfaces for the driven RM become 



The two potentials are displayed in fig. [T] We see that in this ultrastrong coupling regime {g > ^/iij), the 
lower adiabatic potential V^{x) has a double- well structure. This symmetric structure reflects the Z2 parity 
symmetry, which implies that for the double-well potential the spin states are "opposite" between the two 
potential wells. The driving causes the double-well to be asymmetric, and hence the Z2 symmetry is broken. 
The az term in the Hamiltonian opens up a gap between the two potentials (see the inset of the figure). Around 
this avoided crossing, the BOA is likely to break down and it is no longer possible to think about the system 
as two decoupled potentials. For A = 0, the double-well potential is symmetric and for large couplings g the 
spectrum is to a good approximation degenerate for energies E < 0. For positive and moderate energies, this 
quasi degeneracy is lost. These are properties also shared by the Dicke model and there the double- well structure 
characterizes the Dicke phase transition and the corresponding spontaneous breaking of the Z2-symmetry [44| . 
For even larger energies, the anhorminicity deriving from the spin-boson coupling becomes extremely weak 
and the two potentials are approximately harmonic. For a large driving, i.e. X > g, the asymmetry of the 
double- well potential is distinct, which will split the quasi degeneracy. Nevertheless, provided that g is large 
the negative energies can be approximated with those of two harmonic oscillators. Taking all these aspects into 
account, we draw the conclusion that in order to find any non-trivial spectrum the level statistics should be 
explored for moderate and positive energies. This has also been confirmed numerically, i.e. the largest deviation 
from Poissonian statistics is regained in this energy regime. In fig. [2] we show the distribution P{S) of the 
driven RM for energies < -E < 250 and for the same parameters as in fig. [Ij The pronounced "clustering" 
clearly demonstrate the absence of Poissonian statistics. The clustering at small S is even indicating some level 
repulsion. 

As a remark on level statistics. It can be shown that the RM is deeply connected to the E x e Jahn- Teller 
model 14^. While the E x (/3i -I- (32) model shows full blown quantum chaos [4g|, the E x e model displays 
classical chaos and some 'incipience' of quantum chaos [43 . 

6. Level crossings. Parameter dependence of the spectrum of the RM was studied in [48|. In contrast to the 
solvable Jaynes-Cummings model [29|, the energies of the RM show typically avoided crossings within the two 
parity sectors. The driving breaks the Z2 parity and thereby also the crossings arising from this symmetry. 
We have numerically checked this statement, namely that the driving split the crossings between energies with 
different parities. Furthermore, in fig. [2] we already saw some tendencies of level repulsion. Hence, also according 
to this definition, the driven RM seems quantum non-integrable. 

7. Solvability. As we argued above, the question whether the RM is exactly solvable or not is still open. In ref. |34| . 
the conclusions is that the RM is only quasi-exactly solvable. This means that some properties, but not all, are 
obtainable analytically. Note that solvability of the RM does not automatically imply solvability of the driven 
RM. 

Summarizing, according to the standard definitions of quantum integrability the driven RM should not be considered 
integrable. This said, it does not mean that the driven RM is not integrable. As long as there is no accepted definition 
for quantum integrability or exists a solution of the driven RM we simply do not know if it integrable or not. 

B. Thermalization of the Driven Rabi Model 

Frequent in the literature, either non-integrability or chaos are taken as prerequisites for quantum thermalization. 
The previous sections made clear that care needs to be paid for such statements. We have here a model that in many 
aspects should be considered non-integrable and its semi-classical counterpart is chaotic, and still, as we will show, 
we find no evidences for thermalization. 

The numerics is carried out using diagonalization of the truncated Hamiltonian. The truncation in the compu- 
tational basis {|n, 1), |?t., 2)} consists in having an upper limit n < Nn- of the number of bosons. Ntr is taken such 
that our results have converged, i.e. do not depend on A'tr- As local observables we have considered fi = a^d, x, p 
and aa (a = x, y, z), and for non-local observables the "interaction energy" x&x- We will only present statistics of 
the boson number n{t) = {'ijj{t)\n\'ip{t)). Similar results are obtained for the other observables. All our simulation 
of out-of-equilibrium dynamics emerges from a quantum quench. We prepare the system in the ground state of one 
Hamiltonian Hq and at time t = we suddenly shift the parameters of the Hamiltonian to a new one H under 




Figure 3: The boson variance Sn (a) and the scaled boson variance S„/{n}T- The parameters are the same as in fig. [T] 
Non-vanishing variance is a manifestation of non-equilibration. 

which the state evolves. The "initial" Hamiltonian Hq is the RM with g = 0.1 (and thus A = 0), while the system 
Haniiltonian H of eq. ^ typically has g > 1 in order to be in the highly anharmonic regime and A 7^ in order to 
break the Z2 symmetry. The initialized state \ip{t = 0)) is predominantly populating states around the zero energy, 
('0(O)|iJ|'0(O)) ~ 0. In this respect, the eigenstates forming the evolved state are from the irregular part of the 
spectrum in order to maximize the thermalization effect. 




Figure 4: The distribution P{n} for the driven RM. The Gaussian shape signals an incommensurability of the eigenvalues _E^. 
The parameters are the same as for fig. [T] 

Quantum thermalization implies that the boson variance 
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(10) 



vanishes up to order 0{r/^). The coupling dependence of the variance is displayed in fig. [3] (a) for parameters as 
in figs. [T] and [5] For small coupling values g there is some complicated g-dependence, while for larger values the 



variance Sn 



9 



One could imagine that the increased variance for larger g's derives from larger number n{t) 



of bosons. In order to check that this is not the case we show in fig. [3] (b) the scaled variance 5^/ {n)T where 
{n)T = hmT->cx3 y /o dtn{t) is the averaged boson number. Even the scaled variance does not seem to approach zero 
but some finite value for large couplings. In a mean-field approach we can understand why the scaled variance goes 



8 

towards some non-zero value. Within the BOA and deep in the ultrastong couphng regime the ground state of the 
driven RM wih be a coherent state with amplitude a corresponding to the minimum of the lower adiabatic potential 
V~^{x) [43, |4^. For large couplings g > to, \, the coherent amplitude a ^ x ~ g/^f^ so that {n)T ^ g^, and since 
both 5n and (n)T scale as the square of the coupling their ratio should be constant. 
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Figure 5: The Wigner distribution W(x,p) of the field state pi(t) after a time t 

fig-ID 



500 000. The parameters are the same as for 



We continue analyzing the boson statistic by recalling a result by Kac [5^. Given a set of real values {A^} that are 
incommensurate, that is '^^ n^Xi, 7^ for any integers Ui, (except the trivial case n^ = V v), we form the function 

Su{t) — \/^ X]^'=i cos(Aji). The function S,y{t) has a normalized time average S'^{t) = 1. Letting i^ — )■ 00, Kac proved 
that the probability to find Soo{t) between two values a and b is Gaussian, i.e. 



P{a < 5oo(i) <b} = 



dx f 



^x^/2 



(11) 



From this we expect that for incommensurate eigenvalues Ei,, h{t) should be Gaussian. Thus, sampling n{t) at time 
instants {t^} would result in a normal random distribution. For the same initial state as in previous figures, we have 
verified this randomness for the driven RM by calculating the distribution P{n} as shown in fig. 2] The fit to a 
Gaussian with mean {n)T and variance 5„ is almost perfect. Interes ting ly, the Gaussian distribution has also been 
verified for the Jaynes-Cummings model which is definitely integrable 5lJ. Thus, Gaussianity in this respect does not 
prove non-integrability nor chaos. 

One signature for thermalization is that the evolved state \ip{t)) is ergodic and shows seemingly irregular phase- 
space structures [Tj]. For the reduced density operator of the boson field, /5f(t) — ^j=i 2 IOIV'(0)P: "^^ introduce the 
Wigner distribution |5a | 



W{x,p,t)^- I dy{x~y/2\pi{t)\x + y/2)e'Py. 
n 



(12) 



The Wigner distribution is normalized and the marginal distributions agree with the quadrature distributions of 
the boson field. It is not, however, a proper probability distribution since it is not positive definite. One peculiar 
property of the Wigner distribution, also demonstrating that it is not a good probability distribution, is that sub- 
Planck structures are allowed [5j| . In fig. [5] we show an example of the evolved Wigner distribution for the same 
parameters as earlier figures. The time is chosen such that the "collapse" of the initially localized distribution has 
occurred long before the time of the plot. What becomes clear is that the Wigner distribution still shows regular 



interference structures which is expected for non-chaotic time evolution. We have also calculated the corresponding 
Wigner distributions for eigenstates of the driven RM for energies around E rv 0. Surprisingly, these eigenstate 
Wigner distributions are more irregular and they all seem ergodic. Despite this ergodicity of the eigenstates, fig. [5] 
demonstrates that the state does not thermalize. 




Figure 6: The inverse partition ratio rj^ as a function of the coupling g. The remaining parameters are as in fig. [T] 

All numerical results so far suggest that the driven RM does not show quantum thermalization. However, one 
may argue that: (i) only a specific initial state has been considered, and (ii) the driven RM is not a many-body 
model and absence of thermalization could stem from too few contributing states of the sum ([Ij. In order to rule 
out the first possibility, we have checked for several different initial states. In principle, for a system that thermalizes 
the expectations (A) should not depend on details of the initial state but only depend weakly on the system energy 
E. We have thereby focused on analyzing various initial states with different energies. Only states with E > are 
interesting since this is were the spectrum is the most irregular. .E-dependence in (5„ is indeed found, and in all our 
numerical simulations we encounter large fluctuations in n{t). Thus, we can rule out option (i). To get a feeling for 
the finite size effects of our simulation we calculate the inverse partition ratio (|4]) for different couplings g and the 
same type of initial quenched states. The results are shown in fig. |6l As expected, rj^ increases for large couplings. 
If the absence of thermalization derives from finite size effects we should have a decrease of 5n/ {n)T for increasing g 
since corrections from zero should scale as 1/??^. This is not what fig. [3] (b) suggests and we thereby cannot explain 
the large fiuctuations in the variance (5„ as a result of finite size effects. 

IV. CONCLUDING REMARKS 

By considering the driven RM we discussed some ambiguities of quantum integrability and thermalization. Following 
the most commonly used definitions of quantum integrability, the driven RM would be classified as non-integrable. 
The fact that there have been claims that the driven RM is solvable [33i] strengthen the knowledge that quantum 
integrability is a subtle issue. The solvability of the RM, yet alone the driven RM, has however been questioned J34| . 
Instead of being exactly solvable, it is argued that only part of the solutions are analytically obtainable, i.e. the model 
is quasi-exactly solvable. As a non exactly solvable model, a natural exploration is whether the driven RM quantum 
thermalizes. All our numerical simulations indicated that the model do not thermalizes. This, on the other hand, 
proposes that quantum non-integrability is not a necessity for quantum thermalization. Our findings also hint that 
classical chaos cannot be taken as a requirement for quantum thermalization. 

It would be interesting to pursue similar analyses for other models that are in some sense quasi solvable. One 
example would be the Heisenberg XYZ spin-1/2 chain. This model only constitute discrete symmetries, but some 
results, like the ground state energy, can be obtained analytically [5^]. Whether this quasi solvability implies lack 
of quantum thermalization is not known. For the XYZ chain including an external field [5a] there exists no known 
solutions, and thermalization properties of the XY Z model might therby change in the presence of a field. 
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